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Abstract 



We focus on preventing collisions between debris and debris, for which there is no current, effective mitigation strategy. 
We investigate the feasibility of using a medium-powered (5kW) ground-based laser combined with a ground-based 
telescope to prevent collisions between debris objects in low-Earth orbit (LEO). The scheme utilizes photon pressure 
alone as a means to perturb the orbit of a debris object. Applied over multiple engagements, this alters the debris orbit 
sufficiently to reduce the risk of an upcoming conjunction. We employ standard assumptions for atmospheric conditions 
and the resulting beam propagation. Using case studies designed to represent the properties (e.g. area and mass) of 
the current debris population, we show that one could significantly reduce the risk of nearly half of all catastrophic 
collisions involving debris using only one such laser /telescope facility. We speculate on whether this could mitigate the 
debris fragmentation rate such that it falls below the natural debris re-entry rate due to atmospheric drag, and thus 
whether continuous long-term operation could entirely mitigate the Kessler syndrome in LEO, without need for relatively 
expensive active debris removal. 
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1. Introduction 

The threat of catastrophic or debilitating collisions be- 
tween active spacecraft and orbital debris is gaining in- 
creased attention as prescient predictions of population 
evolution are confirmed. Early satellite environment dis- 
tribution models showed the potential for a runaway "Kess- 
ler syndrome" of cascading collisions, where the rate of de- 
bris creation through debris-debris collisions would exceed 
the ambient decay rate and would lead to the formation of 
debris belts ( [Kessler fc Cour-Palais 1978). Recorded col- 
lisions events (including the January 2009 Iridium 33 / Cos- 
mos 2251 collision) and additional environmental modeling 
have reaffirmed the instability in the LEO debris popula- 
tion. The latter has found that the Kessler syndrome is 
probably already in effect in certain orbits, even when the 
models use the extremely conservative assumption of no 



new launches ( Liou & Johnson 2008 , 2009 ) . 



In addition to the UN COPUOS's debris mitigation 
guidelines, collision avoidance (COLA) and active debris 
removal (ADR) have been presented as necessary steps 
to curb the runaway growth of debris in the most con- 
gested orbital regimes such as low-Earth sun synchronous 
orbit (Liou & Johnson, 2009). While active spacecraft 



COLA does provide some reduction in the growth of de- 
bris, alone it is insufficient to offset the debris-debris col- 



lisions growth component (Liou, 2011). Liou & Johnson 



(2009) have suggested that stabilizing the LEO environ- 
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ment at current levels would require the ongoing removal 
of at least 5 large debris objects per year going forward 
(in addition to a 90% implementation of the post mission 
disposal guidelines). Mission concepts for the removal of 
large objects such as rocket bodies traditionally involve 
rendezvous, capture and de-orbit. These missions are in- 
herently complex and to de-orbit debris typically requires 
A'u impulses of order lOOm/sec, making them expensive 
to develop and fly. Additionally, a purely market-based 
program to solve this problem seems unlikely to be forth- 
coming; many satellite owner /operators are primarily con- 
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cerned with the near term risk to their own spacecraft and 
not with long term trends that might endanger their op- 
erating environment, making this a classic "tragedy of the 



Levit & Marshall (2010) provide details of ongoing con- 



commons" (Hardin, 1968j). The cost/benefit trade-off for 
active removal missions makes them unlikely to be pur- 
sued by commercial space operators until the collision risk 
drives insurance premiums sufficiently high to warrant the 
investment. 

To quantify this risk one can look to an example: ESA 
routinely performs detailed conjunction analysis on their 



ERS-2 and Envisat remote sensing satellites (IKlinkrad et 



al. 



2005). Although the number of conjunctions predicted 



annually for Envisat by ESA 's daily bulletins is in the hun- 
dreds, only four events had very high collision probabilities 
(above 1 in 1,000). None of these conjunctions required 
avoidance maneuvers after follow-up tracking campaigns 



reduced orbital covariances, or uncertainties (Klinkrad 



2009). While several maneuvers have been required since 



then, the operational risk is still insufficient to provide in- 
centive for large scale debris remediation effort and this 
highlights the need for low-cost, technologically mature, 
solutions to mitigate the growth of the debris population 
and specifically to mitigate debris-debris collisions which 
owner /operators can not influence with collision avoid- 
ance. Governments remain the key actors needed to pre- 
vent this tragedy of the commons that threatens the use 
of space by all actors. 

Project ORION proposed ablation using ground-based 



lasers to de-orbit debris (fCampbell 1996). This approach 



requires MW-class continuous wave lasers or high energy 
pulses (of order 20 kJ per 40ns pulse) to vaporize the de- 
bris surface material (typically aluminum) and provide 
sufficient recoil to de-orbit the object. ORION showed 
that a 20 kW, 530 nm, 1 Hz, 40 ns pulsed laser and 5 m 
fast slewing telescope was required to impart the Av of 
100-150 m/sec needed to de-orbit debris objects. This was 
technically challenging and prohibitively expensive at that 



time (Phipps et al. 1996). Space-based lasers have also 



been considered, but ground-based laser systems have the 
advantage of greatly simplified operations, maintenance 
and overall system cost. 

In this paper we propose a laser system using only pho- 
ton momentum transfer for debris-debris collision avoid- 
ance. Using photon pressure as propulsion goes back to 
the first detailed technical study of the solar sail concept 
( Garwinj [l£58). The use of lasers to do photon pressure 
propulsion was first proposed by Forward (1962). For the 



application of this to collision avoidance, a Av of Icm/s, 
applied in the anti-velocity direction results in a displace- 
ment of 2.5 km/day for a debris object in LEO. This along 
track velocity is far larger than the typical error growth of 
the known orbits of debris objects. Such small impulses 
can feasibly be imparted only through photon momentum 
transfer, greatly reducing the required power and com- 
plexity of a ground based laser system. Additionally, this 
reduces the potential for the laser system to accidentally 
damage active satellites or to be perceived as a weapon. 



junction analysis research at NASA Ames Research Cen- 
ter, including all-on-all conjunction analysis for the publi- 
cally available U.S Strategic Command (USSTRATCOM) 
Two Line Element (TLE) catalog and simulated future 
catalogs of up to 3 million objects on the Pleiades su- 
percomputer. Their paper also presents early results sug- 
gesting that a high accuracy catalog comparable to the 
USSTRATCOM special perturbations (SP) catalog can be 
generated from the publicly available TLEs; sufficiently ac- 
curate to allow collision avoidance with Av in the sub-cm/s 
range. 



This laser COLA scheme was first proposed in Levit 



Marshall (2010) and it is the purpose of this paper to give 



a more detailed analysis. We focus on assessing the effec- 
tiveness of a laser facility for making orbit modifications. 
The system proposed in this paper uses a 5-10 kW contin- 
uous wave laser mounted on a fast slewing 1.5 m optical 
telescope with adaptive optics and a sodium guide star, 
which allows the laser beam to be continuously focused 
and directed onto the target throughout its pass. 

We start by discussing the underlying physical phenom- 
ena, then describe the baseline system and the design of 
our case study. We conclude by presenting the results of 
a case study, summarizing the potential applications and 
identifying further research. 

2. Methodology: Perturbing LEO debris orbits with 
Radiation Pressure 

In order to assess the feasibility of a collision avoid- 
ance scheme based on laser applied radiation pressure, we 
simulate the resulting orbit perturbations for a number 
of case studies. The laser radiation adds an additional 
force to the equations of motion of the irradiated piece of 
debris, which are then evaluated by a standard high preci- 
sion orbital propagator. Application of a small Av in the 
along- track direction changes the orbit's specific energy, 
thus lowering or raising its semi-major axis and changing 
its period (illustrated in Fig[T]). This allows a debris ob- 
ject to be re-phased in its orbit, allowing rapid along-track 
displacements to grow over time. 

Comprehensive all-on-all conjunction analysis would i- 
dentify potential debris-debris collisions and prioritize them 
according to collision probability and environmental im- 
pact (a function of object mass, material, orbit, etc.), as 
well as screening out conjunctions for which the facility is 
unable to effect significantly (e.g. one involving two very 
massive or two very low A/M debris objects). For con- 
junctions with collision probabilities above a certain "high 
risk" threshold (say 1 in 10,000) we would then have the 
option of choosing the more appropriate object (typically 
lower mass, higher A/M) as the illumination target. Ob- 
jects of lower mass will be perturbed more for a given force 
per unit area. Below we discuss how to approximate the 
area to mass ratio of the object and how to model the 
displacement that is possible with a given system. 
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Figure 1: Schematic of laser system and operations 

2.1. Assessing radiation pressure 

Radiation pressure is a result of the photon momentum. 
If a piece of debris absorbs or reflects incoming photons, 
the momentum transferred leads to a small, but significant 



force. As described in the literature (Mclnnes 1999), the 



resulting force per unit area, i.e. the radiation pressure, is 



F/A = Cry<P=Cry<I/c 



(1) 



where A is the illuminated cross section, / is the intensity 
of the radiation, Cr is the radiation pressure coefficient 
of the object and c is the speed of light. Cr can take 
a value from to 2, where = means the object is 
translucent and Cr = 2 means that all of the photons are 
reflected (i.e. a flat mirror facing the beam). An object 
which absorbs all of the incident photons (i.e. is a black 
body) has Cr = 1. For constant intensities, the resulting 
force can be obtained by simple multiplication. However, 
for larger pieces of debris, the intensity will vary over the 
illuminated cross section. Hence, we choose to implement 
a more accurate description for our simulation, integrating 
over the illuminated cross-section. 



Cr/c j /(x, y) dA 



(2) 



The intensity distribution /(x, y) at the piece of de- 
bris depends on the employed laser, its output power and 
optics, and the atmospheric conditions between the laser 
facility and the targeted piece of debris. In the simplest 
case, I{x^y) will be axisymmetric / = /(r) and follow a 



Gaussian distribution (Siegman, 1986) 



imr) = Ioe 



-2r^ lw{LY 



(3) 



where /q is the maximum intensity of the beam and w is 
the beam width, defined as the radius where the intensity 
drops to 1/e^ of the maximum /q in a given plane at a 
distance L from the laser, /q depends on the beam width, 
as a larger beam width will lead to the energy being dis- 
tributed over a larger area. The beam width is a function 
of the distance L between the laser and the debris, w is 



somewhat controllable but depending on the laser, its op- 
tics and atmospheric conditions, there is a lower limit for 
the beam width. 

The lower limit for it;o(min) foi" an ideal laser propagat- 
ing in a vacuum is given by the diffraction limit. 



^O(mm) ^ \LID 



(4) 



where A is the wavelength of the laser, D is the diameter 
of the focusing optic and L is the distance between the 



optic and the piece of debris (Siegman 1986, p. 676). 

Assuming an object in a 800 km orbit, passing directly 
overhead a station which uses a solid state laser with a 
wavelength of 1pm and a focusing optic with a 1.5 m di- 
ameter, a minimum beam width of 0.6 m would result. 
Increasing the beam width is always possible, but in order 
to maximize the force applied, we assume the beam width 
is at a lower limit. 

In the case of a real laser facility the atmosphere has 
two major effects on beam propagation. First, different 
constituents will absorb and/or scatter a certain amount 
of energy. Second, atmospheric turbulence leads to lo- 
cal changes in the index of refraction, which increases 
the beam width significantly. In addition, the resulting 
time-dependent intensity distributions might not resem- 
ble a Gaussian at all. However, laser engagements in our 
case will take place over time frames of minutes so we 
adopt a time-averaged approach. As common in this field, 
we choose an extended Gaussian model, where the mini- 
mum beam width is increased by a beam propagation fac- 
tor, leading to a reduced maximum intensity. It has been 
shown that this "embedded Gaussian" approach is valid 
for all relevant intensity distributions, allowing simplified 
calculations (Siegman 1991). Even if the Gaussian model 



might not resemble the actual intensity distribution, the 
approach ensures that the incoming time-averaged total 



intensity is correct (Siegman et al. , 1998). The resulting 



intensity at a distance L from the laser depends on the 
conditions on a given path L through the atmosphere. 



/(L,r) =5'snm(^) X r(L) x 



2P 



7rm 



O(min) 



X exp -2Ssum{L) X 



(5) 



O(min) ^ 



where P is the output power of the laser and i^^o(min) is 
the minimum beam diameter in a distance L calculated ac- 
cording to equation |4j This lower limit is increased by the 
Strehl factor Ssum- The total transmitted power is reduced 
by a factor r, accounting for losses through scattering and 
absorption, r and Sgum depend on the atmospheric path 
and this path changes during the engagement as the debris 
crosses the sky. r and Ssum are calculated for each time 
step by integrating atmospheric conditions along the path 
at that time. We use the standard atmospheric physics 
tool MODTRAN 4 (Anderson, 2000) to calculate r. Ssum 
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is a cumulative factor that includes the effects of a less 
than ideal laser system and optics in addition to turbu- 
lence effects. To assess turbulence effects we use the Rytov 
approximation. The Rytov approximation is a statistical 
approach commonly used in atmospheric optics that com- 
bines a statistical turbulence model and perturbation the- 
ory to modify the index of refraction in the wave equation. 
The theoretical background and details of our numerical 
approach are described elsewhere ( |Stupl &: Neuneck 2010 



appendix A), (Stuply^2008^ chapter 2), including additional 
references therein on atmospheric optics and turbulence. 

Our calculations show that turbulence reduces the ef- 
fectiveness of the system by an order of magnitude - prin- 
cipally by increasing the effective divergence. To counter 
those effects, we assume that an adaptive optics system 
with an artificial guide star is used. Such a system mea- 
sures the effects of turbulence and counters them using 
piezoelectric deformable mirrors. The correction has to be 
applied in real time, as local turbulence changes rapidly 
and the guide star moves across the sky as the telescope 
tracks the target. Adaptive optics performance varies de- 
pending on the degree of turbulence in the path of the 
beam and the technical capabilities of the adaptive optics 
system. 

Physical properties of space debris objects vary and for 
a majority of objects some parameters are unknown. This 
makes accurate modeling difficult. A discussion of the key 
parameters and our assumptions follows. 

2.2. Area to Mass ratio 

The acceleration from photon pressure on a debris tar- 
get is proportional to the object's area and inversely pro- 
portional to its mass. To accurately model the photon 
pressure from a beam of width w on an object, both area 
and mass need to be independently known. Since this re- 
search presents an initial feasibility investigation, the di- 
mensions for a random set of debris objects can be in- 
ferred from statistical data on debris size. The ESA Mas- 
ter model provides statistics on observed characteristic size 
distributions (shown in Fig[2| for objects in our region of 
interest, namely sun-synchronous LEO - the most prob- 
lematic region for debris-debris collisional fragmentation 



(Oswald et al., 2006). 



Launch and Mission Related Objects, including rocket 
upper stages and intact satellites, greatly dominate the to- 
tal mass of objects in LEO and are generally too massive to 
be effectively perturbed using photon pressure alone. The 
implication is that this scheme, as presented, would likely 
be ineffective at preventing collisions between two massive 
objects such as rocket bodies or intact spacecraft. How- 
ever, over 80% of all catalogued objects in sun-synchronous 
LEO are debris resulting from explosions or collisions, and 
a significant proportion of these may be effectively per- 
turbed using photon pressure alone since fragments typi- 
cally have high A/M ratios and low masses ( Anselmo fc, 
2010| ). The efficacy of the laser photon pressure 
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Figure 2: MASTER2005 spatial density in sun synchronous orbit 
between 600 km and 1100 km altitude. Note that Launch and Mis- 
sion Related Objects (LMRO) include active, maneuverable satel- 
lites. Additionally, this figure does not include the Fengyun IC and 
Iridium/Cosmos debris (we are awaiting new MASTER data) 



depends on the proportion of collisions that involve high 
area to mass ratio objects in general, and debris fragments 
in particular. The ballistic drag coefficient, defined as the 
product of the dimensionless drag coefficient Cd and the 



area to mass ratio A/M, for an object is given by ( Vallado 



et al., 2007): 



B = CdX A/m = 12.7416215'^ 



(6) 



where (BSTAR) is a free parameter of the orbit deter- 
mination process used to generate TLEs. This relationship 
holds for an atmospheric model that does not vary with 
solar activity but in the case of low solar activity Eq. |6] 
systematically underestimates the ballistic coefficient for 
debris fragments, sometimes by multiple orders of magni- 



tude (Pardini & Anselmo, |2009) . Additionally, the difficul- 
ties in tracking irregular and small debris objects suggests 
that B* for debris objects is less accurate than for large 
rocket bodies or satellites. In fact, a number of objects 
were found in the catalog with no B^ information at all. 

A more accurate method for determining the ballistic 
coefficient is to rescale B by fitting the observed decay of 
the semi-major axis of the object over a long period, us- 
ing an accurate atmospheric model and a high accuracy 
orbit integrator (Pardini & Anselmo 2009|). We imple- 



mented this method by downloading 120 days of TLEs for 
each debris object and then using a standard high preci- 
sion orbit propagator to fit the ballistic coefficient to the 
observed decay of semi-major axis. Assuming Cd = 2.2, 
a reasonable value for the A/M ratio of an object can be 
estimated. 



approach as a long term debris remediation tool therefore 
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2.3. Spin state and Reflectivity 

The spin state of a debris object introduces a degree 
of randomness into calculating the response to directed 
photon pressure. The momentum transferred from ab- 
sorbed photons will be in the incident beam direction. For 
a tumbling target the force vector due to reflection will 
be varying during the engagement, since there will be a 
component of the force orthogonal to the laser incidence 
vector, and for most targets the laser will also induce a 
torque about the center of mass, which we ignore for the 
present. We follow the ORION study and assume that col- 
lision and debris fragments above 600 km will be rapidly 
spinning ( [Phipps et al. 1996). On average, for quickly 
tumbling objects, orthogonal force vectors (due to spec- 
ular reflection) will be zero and the net force vector due 
to diffuse reflection will be directed parallel to the laser 
beam. 

Mulrooney and Matney (2007) suggest that debris has 
global albedo value of 0.13 which in the general case would 
give Cr = 1.13. However, we make a conservative assess- 
ment and neglect the effects of diffuse reflection, assuming 
a force parallel to the laser beam according to equation 
[2| where Cr = 1.0. In reality, the resulting net force will 
likely be larger and for slowly spinning objects the net 
force will not be in the beam direction. 

In an operational setting, one would propagate forward 
a range of laser vector Av (associated with unknowns in 
Cr, A/M etc.) and a range of orthogonal Av to account 
for uncertainties in object forms and spin states. The im- 
plications of the engagement could then be assessed using 
the resulting error ellipsoid of the maneuver e.g. to ensure 
that the maneuver would not cause future conjunctions 
with other objects in the debris fleld. For the purpose of 
this study we also assume that the illuminated cross sec- 
tional area is equal to the effective average cross section, as 
determined by our long term estimation of the drag area. 
This is equivalent to approximating the rapidly tumbling 
object as a sphere of radius equal to this average drag area. 

2.4- Implementation 

For determining the ballistic coefficient of an object 
from the decay of its semi-major axis we used AGI's Satel- 
lite Tool Kit (STK) and an iterative differential corrector 
to flt a high precision orbit to the object's historical TLEs. 

We developed a model for laser propagation in an at- 
mosphere as per Section 2.1 using MATLAB and MOD- 
TRAN 4. Target objects were propagated using a high 
precision propagator in STK, accounting for higher-order 
gravitational terms, a Jacchia-Roberts atmospheric model, 
observed solar flux and spherical solar radiation pressure. 
Laser engagements were modeled by utilizing the MATLAB- 
STK scripting environment, allowing the evaluation of the 
laser intensity and resulting photon pressure at each time 
step. 



3. Baseline System 

Past studies have looked into active debris removal us- 
ing laser ablation. While these favorably assessed the fea- 
sibility of the approach, none of those systems have been 
developed and tested. One reason for this is their reliance 
on what are traditionally military-class systems. These 
are generally not commercially available or are one-of-a- 
kind experimental systems, making them very expensive 
and difficult to obtain. To avoid those shortfalls, we chose 
to restrict this study to medium power commercially avail- 
able lasers and to shorten development times and reduce 
overall cost we also restrict this study to commercially 
available off-the-shelf technology for other parts of the sys- 
tem where possible. Below we outline an example system 
that might be developed today at reasonable cost and the 
following case studies aim to assess whether collision avoid- 
ance is still possible with such a system. 

3.1. Laser 

The intensity that can be delivered to the target (de- 
scribed by equations [4] and [5| is proportional to the laser 
power and inversely proportional to the wavelength. The 
beam quality describes how well the laser beam can be 
focused over long distances, critical for targeting small 
debris objects. Atmospheric transmittance and technical 
constraints puts restrictions on useful wavelengths. For 
targeting sun-synchronous objects the ideal laser facility 
location would be close to the poles and so the equip- 
ment should be low maintenance and ruggedized. Com- 
bining these requirements, and restricting our choice to 
lasers commercially available, we identified an IPG single 
mode fiber laser with a 1.06 pm wavelength. It is electri- 
cally powered with no parts requiring alignment (or that 
can become misaligned) and is designed for 24/7 indus- 
trial applications. The beam quality of this laser is close 
to the diffraction limit (M^ = 1.2) and the output power 
is adjustable up to 5kW (IPG, 2009). 

IPG also manufacturers a 10 kW version and better 
results can be obtained with this higher output power. 
This gives some latitude for the other parameters as dou- 
bling the output power is still possible, albeit at a higher 
cost. As an additional benefit, this low power (compared 
to military systems) makes the system's application as an 
anti-satellite weapon unlikely and thus avoids some of the 
potential negative space security implications. 

3.2. Beam Director and Tracking 

The laser is focused onto the debris using a refiecting 
beam director. The beam director will most likely be an 
astronomical class telescope, potentially modified to man- 
age the thermal effects of continuous laser operation. Ne- 
glecting atmospheric effects the maximum intensity is pro- 
portional to this telescope's aperture. The beam director 
has to be rapidly slewed in order to track the debris and 
the required tracking tolerances become increasingly diffi- 
cult to maintain as diameter and mass increase. Suitable 
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1.5 m telescopes with fast slew capabilities are commer- 
cially available, for example from the company L3, and so 
we choose 1.5 m as our baseline diameter. 

Tracking accuracy for the L3 telescope is of order of 
10~^ arc seconds, which may not be sufficient for track- 



ing small debris in sun synchronous orbit (L3 Communi 



cations Brashear 2008). Hence additional measures have 



to be taken. For laser satellite communications and di- 
rected energy applications, active tracking / closed- loop 
techniques have been developed which are able keep the 
target in the center of the view once it has been acquired 
(e.g. see [Rikerj (2007)). Acquisition is more difficult and 
satellite laser ranging techniques such as beam widening or 
search patterns will be needed to initially find the target. 
It will probably be necessary to use an imaging telescope 
coupled to the beam director to allow simultaneous guide 
star creation, beam illumination and target imaging for 
acquisition and tracking. The Mt. Stromlo facility op- 
erated by Electro Optic Systems (EOS) near Canberra, 
Australia is able to acquire and track debris of 5 cm size 
up to 3000 km range using a 100 W average power pulsed 
laser and a 1.8 m fast slew beam director (Smithl 2007). 



This demonstrates that the target acquisition and tracking 
requirements can be met, although it may prove necessary 
to include a pulsed laser in the proposed system to allow 
for range filtering during target acquisition (as is done by 
SLR systems). 

3.3. Adaptive optics 

Restricting the laser system to a single 5-10 kW facility 
means that sufficient laser intensities can only be reached 
if the effects of atmospheric turbulence are countered by 
adaptive optics. The effectiveness of such a system will 
depend on the turbulence encountered and the technical 
capabilities of the system. 

In our calculations, we assume that the systems ca- 
pabilities for turbulence compensation are comparable to 
the system used in 1998 benchmark experiments (iHiggsl 



1998 Billman et al. 1999), which were conducted to test 



the proposed adaptive optics for the Airborne Laser mis- 
sile defense project. The American Physical Society has 
compiled those results into a relationship of Strehl ratio 
vs. turbulence (Barton et al. 2004 p. 323) and we use 



this relationship in our numerical calculations to set the 
upper limit of the assumed adaptive optics performance. . 
While the ABL is a military system (and has much greater 
output power than necessary for COLA), the Large Binoc- 
ular Telescope has shown a similar performance, reaching 
Strehl ratios up to 0.8( |MPIA[ [2QTq1 ). 

Turbulence effects must be measured in order to be 
compensated.. We assume that a laser guide star (po- 
sitioned ahead of the target to account for light travel 
times) is used as a reference point source and compen- 
sation for de-focus and higher order turbulence effects is 
ideal. However, tip/tilt correction requires a signal from 
the real object and not the guide star. We calculate the 



negative impact of this so called tilt-anisoplanatism and 
lower the intensities accordingly. 

3.4- Location and Atmospheric conditions 

The described system is designed to illuminate debris 
in sun-synchronous orbits, so to maximize engagement op- 
portunities we favor a location as close as possible to the 
poles. Additionally, situating the facility at high altitude 
reduces the atmospheric beam losses and turbulence ef- 
fects. An ideal site would be the PLATeau Observatory 
(PLATO) at Dome A in Antarctica, which is at 4 km alti- 
tude and is in the driest region of the world. For compar- 
ison we also considered Maui and Mt. Stromlo, since they 
already have facilities that might be upgraded to test this 
concept, and a hilltop near Fairbanks, Alaska due to its 
high latitude and ease of access compared to arctic terri- 
tory. 

Atmospheric conditions will have a major impact on 
the performance of the system. Site selection and dome 
design will have to take this into account to minimize losses 
and down time. For this study we chose standard condi- 
tions for turbulence and atmospheric composition (Hufnagel- 
Valley 5/7 turbulence and the US Standard Atmosphere 
(1976), MODTRAN set to 365 ppm CO2, Spring/Summer 
conditions, and 23km surface meteorological range). 

3.5. Scalability 

While the laser parameters are readily available using 
a datasheet, tracking accuracy and adaptive optics perfor- 
mance are less certain. Since the effect of laser engage- 
ments is cumulative, one could both increase the power of 
the laser and use multiple stations, engaging debris from 
different locations, if adaptive optics performance or accu- 
rate tracking becomes more difficult than expected (or if 
one wants to do collision avoidance for lower A/M or heav- 
ier debris objects). For example, by upgrading the laser 
to a 10 kW model and having 3 or 4 facilities the effect of 
this system can be increased by an order of magnitude. 

3.6. Operational Considerations 

In general we want to lower the orbits of debris ob- 
jects to reduce their lifetime so the optimal tasking of the 
laser-target engagement is to begin illuminating the tar- 
get from the horizon and to cease the engagement when 
the target reaches its maximum elevation (simulated en- 
gagements start at 10° elevation to approximate acquisi- 
tion delays). The main components of the net force for 
an overhead pass are in the anti-velocity and radial direc- 
tions. Engaging during the full pass would result in a net 
radial Av, which results in less rapid displacement over 
time from the original trajectory. 

Target acquisition and tracking at the start of each en- 
gagement will produce track data and, if a pulsed laser is 
used for acquisition, ranging data similar to that produced 
by the EOS Space Debris Tracking System ( |Sniith[ |20Q7| . 
This would allow orbit determination algorithms to reduce 
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Figure 3: Orbital decay of semi-major axis for Akari lens cap. The 
"Fitted Orbit" represents the orbit decay using the rescaled A/M 
ratio, as fitted to the TLEs with a highly accurate special perturba- 
tions propagator. 



the error covariance associated witii tiiat object's orbit - 
iielpful for space situational awareness (SSA) in addition to 
down range target re-acquisition. Laser campaigns would 
only need to continue until the collision risk has been re- 
duced to an acceptable level - which can be either through 
improved covariance information and/or through actual 
orbit modification. 



4. Resulting Capabilities 

To quantify the effectiveness of this laser scheme on de- 
bris objects we start by demonstrating our method for an 
object of known mass and area. A discarded lens cap from 
the Japanese Akari IR space telescope was chosen as the 
demonstration object (U.S. Catalog ID: 29054). We nomi- 
nally chose 01 January 2011 00:00:00 UTC as the starting 
time for all simulations. The lens cap is approximately 
a fiattened hemispherical dome of mass 5 kg, with a di- 
ameter of 80 cm and a thickness of approximately 10 cm. 
These parameters represent a large debris fragment. This 
lens cap orbits in a near circular orbit at about 700 km 
altitude, with an inclination of 98.26°. 

Fitting the observed orbital decay of the lens cap over 
120 days (shown in Fig[3| to derive the ballistic coefficient 
gave A/M = 0.04. This is close to the minimum ballis- 
tic coefficient possible with the known object dimensions, 
suggesting that the lens cap has stabilized to present a 
minimum cross-section and to minimize drag forces. We 
initially use this area for radiation pressure calculations, 
even though the surface visible to the laser is likely to be 
larger. 

Fig [4] shows how the beam radius varies due to the 
changing beam path as the lens cap passes over the facil- 
ity, with the engagement ending at the maximum eleva- 
tion. The peak intensity (at the center of the beam) and 



resultant power on the target are a minimum at the low- 
est elevation and increase throughout the 5 minute pass. 
The resulting displacements from 5 kW laser engagements 
during the first half of each pass of the debris object over 
the laser during a 48 hour period (25 engagements in the 
case of PLATO) are compared in Fig [s] for four separate 
locations. The in-track rate of displacement, or velocity 
difference, resulting from the illumination campaign is 82 
m/day. Using the approach given in Levit & Marshall 
( |2010[ ), an analysis of 70 days of TLEs for the lens cap 
showed that the orbital in-track error grows by an aver- 
age of 178 m/day, when propagated with a fitted numeri- 
cal orbit propagator. . This method alone would not be 
sufficient to detect a maneuver on this object. However, a 
10 kW facility would generate 161 m/day which may well 
be detectable. 

For a conjunction of two objects with similar magni- 
tude error to the Akari lens cap (and provided that one 
arranges the engagement geometry so as to increase the 
current predicted miss distance (e.g. by appropriately 
choosing between velocity vector and anti-velocity vector 
nudging)) such a system may be sufficient to significantly 
reduce the collision probability of a conjunction. With 
higher accuracy data based on any of (a) access to the U.S. 
Strategic Command unclassified SP catalog, (b) improved 
orbits obtained from tasked radar /optical tracking or (c) 
TLE improvement scheme proposed in Levit & Marshall 
(2010), it is highly likely that the laser can provide more 
than sufficient Av to overwhelm the orbit/propagation er- 
rors, at least for objects of sufficiently high area to mass 
ratio. As an initial guide point, we will hereafter con- 
sider displacements of more than 200 m/day as significant 
in that they are likely to overwhelm orbit errors associated 
with propagating high accuracy debris orbits. 

We chose a random subset of 100 debris objects from 
the U.S. TLE catalog with inclinations between 97 and 102 
degrees and orbit altitudes between 600 and 1100 km. Our 
selection was limited to this number by the computational 
requirements of running these simulations. Characteristic 
sizes were assigned to these objects to give a representative 
size distribution, shown in Fig[6j 

The A/M ratio of each object was determined (see Fig 
[t]) by rescaling the ballistic coefficient, allowing us to derive 
mass values for the set. 

The mean A/M after rescaling was 0.24m^/kg and the 
median was O.llm^/kg. Two days was selected as a rea- 
sonable minimum conjunction warning lead time, during 
which the laser system could be employed. The laser was 
tasked with illuminating the target for the first half of each 
pass for 48 hours and the resultant displacement (from the 
unperturbed orbital position) was generated for the next 
five days. 

As the size of the object increases beyond the beam 
width, the force on the object asymptotically approaches 
Fmax = CrX 1/c X I^ax X TT X (1/2) xw^ff. There is there- 
fore an upper limit on the mass of an object that can be 
sufficiently perturbed using laser applied photon pressure 
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Figure 4: The behavior of the beam as it tracks the Akari lens cap through a single near-overhead pass. 




Figure 5: Displacement of Akari lens cap from unperturbed orbit after 2 days of laser engagements, plotted for different system locations (for 
details see table [T]). 
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Figure 6: Size distribution for 100 debris objects in sun-synchronous 
LEO, generated using MASTER2005's characteristic size distribu- 
tions. 



16 
14 
12 
10 



o 
O 















































1 








































































































































n 


























mn 





-4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 
LogioA/IVI [m2/kg] 

□A/M rescaled by fitting orbital decay DA / M derived directly from B* 

Figure 7: Debris subset A/M distribution, as inferred by a long term 
(120 day) statistical orbital decay assuming = 2.2. 



with any given system. This hmit depends strongly on 
the geometry of the laser-target interaction, so we do not 
derive this limit analytically. To give an idea of this upper 
mass threshold, objects with masses greater than 100 kg 
were all perturbed by less than lOOm/day. As expected, 
photon pressure is generally not sufficient for maneuvering 
massive objects. 

For a single 5kW laser facility located at PLATO in 
Antarctica, the displacement from the unperturbed orbit 
for 100 objects is plotted in Fig [Sj After a two day laser 
campaign it was found that 43 of 100 objects were di- 
verging from their unperturbed orbit by more than 200 
meters per day and 13 by more than 500 meters per day. 
For a 10 kW laser, 56 objects where perturbed more than 
200 m and 34 more than 500 m. A number of other "suc- 
cess rates" , defined as the number of objects displaced by 
more than x m/day, are shown in Table [l] Situating such 
a laser system in Antarctica may prove infeasible, so for 
comparison the simulation was run for the case of a single 
laser situated at the Air Force Maui Optical and Super- 
computing site in Hawaii, at Mt. Stromlo in Australia 
and at a fictional location near Fairbanks, Alaska. Table 
[l] shows the success rate of the system at these different 
locations for a 5kW and 10 kW laser system. 

Since the targets are all approximately sun synchronous 
the effectiveness of sites away from the polar region is 
greatly reduced, as expected. Mt. Stromlo and Maui show 
similar levels of performance. The additional atmospheric 
losses at Mt. Stromlo's lower altitude are offset by its 
higher latitude. Alaska performs better due to its higher 
latitude, but would benefit from being situated at higher 
altitude. The success rates shown in Table [T] are meant 
to give a qualitative estimate of the campaign's effective- 
ness at avoiding collisions. The true effectiveness of a laser 
campaign is measured by re-evaluating the collision proba- 
bility to determine whether it has decreased sufficiently to 
be confident of a miss. The collision probability is derived 
from the orbital covariance of the two objects, which was 
not available for this analysis. Therefore we do not per- 
form a thorough collision probability analysis, but rather 
present the range displacements resulting from the simu- 
lated laser illumination campaign. 

A 200 m/day range displacement is equivalent to a Av 
impulse of about 0.08 cm/s in the anti-velocity direction. 
Typical Envisat collision avoidance maneuvers have been 
of the order of a few cm/s, but were usually performed 
within a few hours of the conjunction epoch. Satellite 
operators want to minimize a maneuver's impact to the 
lifetime and mission schedule and therefore take the deci- 
sion at the latest possible time to be sure that the maneu- 
ver is actually necessary. Additionally, for remote sens- 
ing satellites where lighting angles are important, maneu- 
vers are often selected to quickly raise or lower the orbit 
to increase the radial miss distance, rather than rephras- 
ing the satellite in True Anomaly, and/or they are com- 
bined with station-keeping maneuvers. For debris-debris 
collision avoidance using a laser this is not a concern and 
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Table 1: Success Rates for 5 and 10 kW laser systems, compared for different sites. The Success Rates are defined as the number of objects 
displaced more than 50, 100, 200, 500 or lOOOm/day. 



Site Parameters Success Rates (daily displacement) 



Power / Location 


Latitude 


Longitude 


Altitude 


50 m 


100 m 


200 m 


500 m 


1000 m 


5kW PLATO, Antarctica 


-80.37 


77.35 


4.09 km 


74 


56 


43 


13 


5 


5kW AMOS, Hawaii 


20.71 


-156.26 


3.00 km 


30 


13 


5 


4 


2 


5kW Mt. Stromlo, Australia 


-35.32 


149.01 


0.77 km 


11 


4 


4 


3 





5kW Eielson AFB, Alaska 


64.85 


-148.46 


0.50 km 


31 


12 


5 


4 


2 


10 kW PLATO, Antarctica 


-80.37 


77.35 


4.09 km 


89 


74 


56 


34 


13 


10 kW AMOS, Hawaii 


20.71 


-156.26 


3.00 km 


42 


30 


13 


5 


4 


10 kW Mt. Stromlo, Australia 


-35.32 


149.01 


0.77 km 


29 


12 


4 


4 


3 


10 kW Eielson AFB, Alaska 


64.85 


-148.46 


0.50 km 


48 


31 


12 


4 


4 
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engagement campaigns may begin much earlier (i.e. two 
days before), letting small changes to the semi-major axis 
re-phase the target over longer periods. 



[Levit fc Marshall ( 2010, ) suggest that batch least-squares 
fitting techniques can generate high accuracy orbital state 
vectors with errors that grow at about lOOm/day. This 
error growth is of the same level as that provided by the 
high accuracy special perturbations catalog(s) maintained 
by the US Strategic Command (Boer s et al.[|2000" ). Given 
either of these sources, a range displacement of 200 m/day 
would dominate the growth of the object's error ellipse and 
would thus likely be sufficient for collision avoidance, but a 
full collision probability analysis is needed to confirm this. 
Additionally, data from initial engagements could reduce 
the size of the error ellipse, meaning that less range dis- 
placement (or, equivalently, less Av) will be required to 
reduce the collision probability. 

5. Discussion on Next Steps and Implications 

5.1. Further Research 

Immediate follow up work should focus on reducing 
the uncertainty of modeling assumptions to improve the 
statistics presented here. Near-term improvements should 
include the following: 

1. Test the effect of this scheme in long-term evolu- 
tionary models, such as the NASA LEGEND model 
( Liou et al.[ 2004|. By considering the long term 



consequences of shielding the "high impact" popu- 
lation (objects of both high collision cross-section 
and large mass) from the type of objects for which 
photon pressure is effective we could determine how 
many objects would need to be shielded to halt the 
cascading growth of debris in low Earth orbit. This 
would provide a better prediction of the long term 
effectiveness of the system. 

2. Radar Cross Section (RCS) data might be used to 
determine the characteristic sizes for individual de- 
bris objects, instead of - as we have done - using 
randomly assigned sizes that match the observed dis- 
tribution. This would allow simulations using more 
accurate object areas and masses. There is some un- 
certainty in the accuracy of RCS measurements, and 
further research and analysis should be conducted 
before adopting this approach. 

3. The simulations should be run for a much larger set 
of objects, and in a wider range of orbital regimes, to 
allow useful statistics to be generated and a metric 
devised to identify the class of objects for which the 
system is truly effective. 

4. Error covariances should be generated for each simu- 
lated object's orbit. This would allow us to estimate 
the change in collision probability resulting from con- 
secutive engagements, a far more useful measure of 
the systems capability than the simple range dis- 
placement. 



5. A systematic parameter optimization study needs to 
be done to identify the best combination of laser 
power versus number of facilities, the optimal loca- 
tions for these facilities, the most advantageous en- 
gagement strategy and the ideal combination of laser 
and optics. 

6. Spin assessment. Since the illuminations can provide 
torques to the debris objects being illuminated, it is 
prudent to research, in detail, the effects that this 
could have. For example, changing the spin rates 
could alter the drag coefficient and make it harder 
to predict the orbit position of that debris object. It 
would also effect the decay lifetime, potentially mak- 
ing it longer. Finally, it could reduce the object's 
radar cross-section. None of these issues seem on 
first analysis to be a significant challenge to the sys- 
tem's overall utility but they demand detailed con- 
sideration. 

7. Finally, the policy implications need consideration. 
These include the problems of debris ownership (and 
potential need for transfer of that ownership) and 
associated liability of maneuvering a piece of debris. 
There are also potential security concerns for the sys- 
tem which may demand solutions similar to laser de- 
confiiction, as practiced by the ILRS (Pearlman et 



al. 2002). 



5.2. Technology Demonstration 

Following the aforementioned further research and a 
comprehensive engineering and costing analysis, a techni- 
cal demonstration would be the logical next step. This 
could most easily be accomplished by integrating a con- 
tinuous wave fiber laser (and adaptive optics if necessary) 
into an existing fast slewing optical telescope and demon- 
strating the acquisition, tracking and orbit modification 
of a known piece of debris (a US-owned rocket shroud for 
example). The thermal, mechanical and optical implica- 
tions of continuous 5kW IR laser operations would need 
to be addressed via engineering simulation first, and prob- 
ably verified in actual tests. Eventual candidates for a 
demonstration include the EOS Mt. Stromlo facility and 
the Advanced Electro-Optical System at AMOS. AEOS 
has demonstrated large-aperture debris tracking with the 
180W HI-CLASS ladar system ( |Kovacs et aL|[2mT| . EOS 
is routinely performing laser tracking of LEO debris ob- 
jects smaller than 10 cm in size from this facility (Greene^ 
20021 . The EOS facility would probably require the fewest 
modifications to incorporate a higher power CW fiber laser 
for a technology demonstration. Since the 5 kW laser costs 
$0.8M, we speculate that the direct cost of adapting such 
a system would be of order $1-2M. In addition, it may be 
possible to perform a near-zero cost demonstration using 
existing capabilities such as those of the Starfire Optical 
Range at Kirtland AFB. It should be noted that the au- 
thors know of no relevant system that already has adaptive 
optics capable of fast slew compensated beam delivery to 
LEO. 
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Having demonstrated the method on an actual piece of 
debris, a fuhy operational system could be designed and 
located at an optimal site, or appended to a suitable ex- 
isting facility. Preliminary discussions with manufacturers 
suggest that the capital cost of the laser and primary beam 
director would be around $3-6M. The cost of the necessary 
primary adaptive optics and tracking systems (including 
secondary lasers and tracking optics) are less clear at this 
stage since there are a number of ways that a working so- 
lution could be engineered. Further engineering analysis 
is necessary before accurate overall system costs can be 
estimated. There is advantage to making the system an 
international collaboration in order to share cost, to ease 
certain legal obstacles to engaging space objects with var- 
ied ownership and to reduce the likelihood of the facility 
being viewed negatively from a security stand point. This 
system would coincidentally complete many of the steps 
(both technical and political) necessary to implement an 
ORION-class laser system to de-orbit debris, potentially 



clearing LEO of small debris in just a few years (Phipps et 



al. 1996), if it was deemed useful to do that in addition. 



A key component for the proposal herein would also be an 
operational all-on-all conjunction analysis system, the cost 
of which is also uncertain but likely to be small compared 
to the other system costs to operate (and which would also 
benefit from including multiple international datasets). 

5.3. Potential Implications for the Kessler Syndrome 



Liou & Johnson ( 2009 ) have identified the type of "high 



impact" large mass, large area objects that will drive the 
growth of the LEO debris population from their catas- 
trophic collisions. In the LEO sun synchronous region 
the high impact debris mass is approximately evenly di- 
vided between large spacecraft and upper rocket bodies 
(Liou 2011). ESA routinely monitors all conjunctions 



with objects predicted to pass through a threat volume 
of 10 km X 25 km X 10 km around its Envisat, ERS-2 and 
Cryosat-2 satellites using their Collision Risk Assessment 
tool (CRASS). These satellites are operational and maneu- 
verable, but their orbit and mass and area profiles' make 
them analogous to Liou and Johnson's high impact ob- 
jects. We therefore use these satellites as a proxy for the 
high impact population. 

75% of conjunctions with Envisat 's threat volume in- 
volve debris (i.e. not mission related objects, rocket bodies 
or other active spacecraft). Significantly, 61% of all En- 
visat conjunctions involve debris resulting directly from 
either the Fengyun 1-C ASAT test or from the Iridium 
33/Cosmos 2251 collision. For ERS-2 and Cryosat-2 (at 



a lower altitude) these figures are similar (Flohrer et al. 



|2009 ) . It is clear that debris resulting primarily from colli- 
sion and explosion fragments is most likely to be involved 
in collisions with large objects in the LEO polar region. 

These statistics suggest that it may be possible to shield 
high impact objects from a significant proportion of catas- 
trophic collisions with less massive debris such as frag- 
ments by using a ground based medium power laser. If 



75% of conjunctions with high impact objects involve de- 
bris (as suggested by Envisat conjunctions) and our anal- 
ysis of 100 random debris objects suggest that 43% can be 
significantly (>200m/day) perturbed using our baseline 
5kW system, then it may be possible to prevent a third 
of all conjunctions involving the high impact population. 
Increasing the laser power to 10 kW would raise this figure 
to 42%. 

Additionally, LEGEND simulations have shown that 
catastrophic collisions involving intacts (spacecraft and rocket 
bodies) and fragments are slightly more likely than colli- 
sions involving only intacts (Liou |2011[ ). Using these col- 
lision statistics, and assuming 200m/day is sufficient to 
insure a clear miss, we see that a single 5 kW system could 
prevent nearly half of all catastrophic collisions involving 
debris fragments, and about 28% of all collisions, includ- 
ing intact-intact collisions. Obviously an intact-intact col- 
lision is a bigger debris source than an intact-fragment of 
fragment-fragment since it involves two massive objects. 
Further LEGEND modeling would be able to quantify the 
degree to which the scheme reduces debris sources. 

Of course one is not limited to shielding one object. We 
posit that it may be possible to use laser photon pressure 
as a substitute for active debris removal, provided a suf- 
ficient number of high impact objects can be continually 
shielded to make the two approaches statistically similar. 
Indeed, the routine active removal of 5 large debris ob- 
jects per year is predicted to prevent 4 intact-intact, and 
5 intact-fragment catastrophic collisions over the next 200 
years (?). With an effective all-on-all conjunction analy- 
sis system to prioritize engagements and considering that 
every engagement reduces the target's orbital covariance 
(thereby halting unnecessary engagement campaigns) it is 
plausible that far more objects may be shielded than are 
required to make the two approaches equivalent in terms of 
preventing the number of catastrophic collisions (a LEG- 
END simulation may confirm this). 

For a facility on the Antarctic plateau the laser would 
be tasked to an individual object for an average of 103 
minutes per day. The laser can only track one target at 
a time, but average pass times suggest that it is possi- 
ble to optimize a facility to engage '^10 objects per day. 
The Envisat conjunction analysis statistics suggest around 
10 high risk (above 1:10,000) events per high impact ob- 
ject, per year ( [Flohrer et al. 2009). If improved accuracy 
catalogs or tracking data become available then it is feasi- 
ble that the system could engage thousands of (non-high 
impact) objects per year, or conversely that up to hun- 
dreds of high impact objects could be shielded by one fa- 
cility per year. This is an order of magnitude more objects 
than one needs to remove in order to stabilize the growth 



(Liou fc Johnson^ 2009). Preventing collisions on such a 
large scale would therefore likely reduce the rate of debris 
generation such that the rate of debris reentry dominates 
and the Kessler syndrome is reversed at low enough alti- 
tudes. Continued operation over a period similar to the 
decay timescale from the orbital regions in question (typ- 
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ically decades) could thus reverse the problem. Addition- 
ally, scaling such a system (eg. multiple facilities) on the 
ground would be low cost (relative to space missions) and 
can be done with currently mature technology, making it 
a good near term solution. Further, if the current analysis 
proves optimistic, raising the power to 10 kW and having 
3-4 such facilities would increase the number of conjunc- 
tions that it is possible to mitigate by a further order of 
magnitude, and also would raise the maximum mass and 
reduce the minimum A/M threshold for the system. 

5. J1.. Additional Applications 

The described system has a number of alternative uses, 
which may further improve the value proposition. 

Firstly, orbit tracks are a byproduct of target acquisi- 
tion that can be used for orbit determination. Correlat- 
ing these tracks would allow the generation of a very high 
accuracy catalog, similar to that being produced by the 
EOS facility at Mt. Stromlo. The return signal from laser 
illumination will potentially provide data for accurate es- 
timation of debris albedo and, if the object is large enough 
to be resolved, size, attitude and spin state; thus helping 
space situation awareness more generally. 

Secondly, the concept of shielding high impact debris 
objects can be applied to protecting active satellites. The 
laser system could begin engaging the debris object follow- 
ing a high risk debris-satellite conjunction alert. The ini- 
tial engagements would provide additional orbit informa- 
tion that may reduce the risk to an acceptable level. Con- 
tinued engagement would perturb the debris orbit, poten- 
tially saving propellant by avoiding the need for a satellite 
maneuver. This could even be provided as a commercial 
service to satellite operators wishing to extend operation 
lifetimes by saving propellant. 

Lastly the laser system may also prove useful for mak- 
ing small propellant-less maneuvers of satellites, including 
those without propulsion, provided the satellite is suffi- 
ciently thermally protected to endure 5-minute periods of 
illumination with a few times the solar constant. This 
could be used to, for example, enable format ion- flying clus- 
ters of small satellites, or perform small station- keeping 
maneuvers. Being able to extend smallsat lifetimes with- 
out launching to higher altitudes or being able to gradually 
re-phase a satellite in True Anomaly may also have com- 
mercial applications. 

6. Conclusion 

It is clear that the actual implementation of a laser 
debris-debris collision avoidance system requires further 
study. Assumptions regarding the debris objects proper- 
ties need refinement and a detailed engineering analysis is 
necessary before a technology demonstration can be con- 
sidered. However, this early stage feasibility analysis sug- 
gests that a near-polar facility with a 5 kW laser directed 
through a 1.5 m fast slewing telescope with adaptive optics 



can provide sufficient photon pressure on many low-Earth 
sun-synchronous debris fragments to substantially perturb 
their orbits over a few days. Additionally, the target ac- 
quisition and tracking process provides data to reduce the 
uncertainties of predicted conjunctions. The laser need 
only engage a given target until the risk has been reduced 
to an acceptable level through a combination of reduced or- 
bital covariance and actual photon pressure perturbations. 
Our simulation results suggest that such a system would 
be able to prevent a significant proportion of debris-debris 
conjunctions. 

Simulation of the long term effect of the system on 
the debris population is necessary to confirm our suspi- 
cion that it can effectively reverse the Kessler syndrome 
at a lower cost relative to active debris removal (although 
quite complementary to it). The scheme requires launch- 
ing nothing into space - except photons - and requires no 
on-orbit interaction - except photon pressure. It is thus 
less likely to create additional debris risk in comparison 
to most debris removal schemes. Eventually the concept 
may lead to an operational international system for shield- 
ing satellites and large debris objects from a majority of 
collisions as well as providing high accuracy debris track- 
ing data and propellant-less station keeping for smallsats. 
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